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Abstract: The "Cold Spot" in the CMB sky could be due to the presence of an anomalous huge 
spherical underdense region - a "Void" - of a few hundreds Mpc//i radius. Such a structure would 
have an impact on the CMB two-point (power spectrum) and three-point (bispectrum) correlation 
functions not only at low I, but also at high i through Lensing, which is a unique signature of a 
Void. Modeling such an underdensity with an LTB metric, we show that for the power spectrum the 
effect should be visible already in the WMAP data only if the Void radius is at least 1/ > 1 Gpc//i, 
while it will be visible by the Planck satellite if 1/ > 500 Mpc//i. We also speculate that this could 
be linked to the high-^ detection of an hemispherical power asymmetry in the sky. Moreover, there 
should be non-zero correlations in the non-diagonal two-point function. For the bispectrum, the 
effect becomes important for squeezed triangles with two very high ts: this signal can be detected 
by Planck if the Void radius is at least L > 300 Mpc//i, while higher resolution experiments should 
be able to probe the entire parameter space. We have also estimated the contamination of the 
primordial non-Gaussianity /nl due to this signal, which turns out to be negligible. 

PACS numbers: 98.80.Cq,98.80.Es, 98.65.Dx, 98.62.Sb 



I. INTRODUCTION 

One of the anomalies present in the recent WMAP [l[ measurement of the Cosmic Microwave Background 
(CMB) is the so-called Cold Spot 0, 1^: a large spherical region on an angular scale of about 10° that 
appears to be anomalously cold. The probability that such a pattern would derive from Gaussian primordial 
fluctuations is estimated to be about 1.8%, using the most conservative criteria 0, While this could 
still be a statistical fluke, some authors have put forward the idea that it could be due instead to some 
object on the line-of-sight between us and the Last Scattering Surface (LSS) 0, 

In a previous paper Q we have explored some observational consequences of the proposed idea that the 
Cold Spot is due to a "Void" 0, Q i e. to an anomalously large underdense region of some unknown origin. 
Traveling through a Void, photons are redshifted due to the fact that the gravitational potential is not 
exactly constant in time, the so-called Rees-Sciama (RS) effect 0- We have thus computed observational 
quantities associated to the RS effect, focusing in particular on the statistical properties of the CMB: the 
two-point (power spectrum) and three-point (bispectrum) correlation functions. 

In the present paper we extend the analysis of by considering also the Lensing effect, namely the 
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deflection that occurs to a CMB photon travchng through a Void. The interesting fact is that this deflection 
is only present if there is something on the hne of sight, while it would be absent if the Spot is just 
a statistical fluke of the primordial large-scale fluctuations, thus representing a unique signature on the 
CMB maps, which can rule out or confirm the existence of such a Void. Moreover the Lensing effect is 
correlated with the RS effect, which can be seen in the some observables (such as the three-point correlation 
function 0]). 

The paper is organized as follows. In section II we consider the lensing produced by an underdense 
region and compute its profile and its decomposition in spherical harmonics. In section III we compute 
the contribution to the two-point function. In section IV we compute the contribution to the bispcctrum, 
and evaluate the signal-to-noise ratio and the possible contamination on a primordial non-Gaussian signal. 
In section V we draw our conclusions. In appendix |X] we discuss how our result change upon inclusion 
of a cosmological constant, while in appendix |B] we review the flat-sky approximation for the primordial 
bispectrum. 

While this paper was in preparation, a similar approach has been taken by who have studied the 
lensing effect from a Void (and also from a texture), showing that high-resolution experiments, which focus 
on a small area of the sky, should be able to detect it. We briefly comment on this result in section III. 



II. A VOID IN THE LINE OF SIGHT: REES-SCIAMA AND LENSING EFFECTS 

As in Q, we consider the following cosmological configuration: an observer looks at the CMB through 
a spherical Void of radius L, located at comoving distance D from us in the direction of the z axis. We 
assume that the Void does not intersect the LSS and that we are not inside it. The observer receives 
from the LSS the CMB photons, whose fiuctuations we assume to be adiabatic, nearly-scale invariant and 
Gaussian - as those generated e.g. by the usual inflationary mechanism. We will also assume that the 
location of the Void in the sky is not correlated at all with the primordial temperature fluctuations coming 
from inflation; this is true, for example, if such structure comes from a different process, such as nuclcation 
of bubbles. For simplicity we disregard the effect of a cosmological constant. We show in appendix 1X1 how 
the expressions can change, after inclusion of a cosmological constant. 

We model the Void's inhomogcneous region via a spherically symmetric Lcmaitrc-Tolman-Bondi (LTB) 
metric, matched to a Friedmann-Lemaitre- Robertson- Walker (FLRW) flat model Our density proflle 
turns out to be the one of a "compensated" Void, i.e. the miderdense central region is surrounded by a 
thinner external ovcrdcnsc region. This follows from the matching conditions, requiring that the Void does 
not distort the FLRW metric at large distance^. Note that, due to the compensation, there will be no 
lensing effect on the photons which travel outside the LTB region. The Void is then characterized by the 
following quantities: i) the distance D from us and its centre, ii) its radius L (namely the radius of the 
inhomogcneous LTB region), iii) the shape of the Void density proflle and iv) the amplitude of the density 
contrast, which we parameterize the value at the centre (which is defined to be negative). Note that, 
since the subtended angle 29 l is small, we have 9l~ L/D. 

Given our configuration, the observer detects one particular realization of the primordial Gaussian per- 



^ For example, if the Void comes from a primordial bubble of true vacuum, this is a physically consistent requirement: in 
fact a bubble would have a thin wall with localized gradient energy, that compensates for the lower energy contained in 
the true vacuum in the interior region. 
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turbations on the LSS plus the secondary effects due to tfiis anomalous structure. As in [6|, we want to 
give the theoretical prediction for the two-point and three-point correlation functions, in order to compare 
with the observations. Wc write the observed temperature fluctuation as the sum of three components: 

Ar(n) _ Ar(n) '-^^ AT{n) '^^^^ ATjn) 
where (P) stands for primordial, (RS) for Rees-Sciama and (L) for Lensing. Each fluctuation is defined as: 



T T 



P.RS.L, (2) 



with the bar representing the angular average over the sky and T = = 2.73A'. The RS temperature 

fluctuation is smaller than the primordial one, but larger than the temperature fluctuation induced by 
Lensing. 



A. Rees-Sciama Temperature Profile 

We now briefly review Q how to compute the shape for AT^^^^ jT . In 0] we have modeled the Void in 
the line of sight as a perturbation of a FLRW metric with a Newtonian potential 3>: 

ds^ = a^{T) [-(1 + 2$)dT2 + (1 - 2^)dx'dx^] , (3) 

where r is the conformal time and where are dimensionless comoving coordinates. The potential $(r) 
is given as a function of the dimensionless comoving radial coordinate r by: 

^i-) = £ mrdf , r<r,^ ^(^^^o , (4) 

while $(r) vanishes for r > r^. Here Hq is the present Hubble constant Hq = h/3000 Mpc~^. The 
arbitrary function k{r) represents the local curvature and determines the shape of the density proflle. 
This approximation is valid as long as k(r) is small. The only constraints on this function come from 
the smoothness of the density profile at the centre (which dictates fc'(O) = 0, where the prime denotes a 
derivative with respect to the r coordinate) and from the requirement that the LTB patch matches to a 
flat FLRW universe {k'{L) = k{L) = 0). We have chosen arbitrarily the function k{r) as follows: 



k{r) = ko 



rL 



(5) 



In the rest of this paper we focus on the case a = 4 but, as shown in [6|, the results for other cases, e.g. 
a = 2, are not qualitatively different. 

The RS temperature fluctuation is effectively described by two parameters: its amplitude at the centre 
of the Void, A = AT{z)'-^^^ /T, and its angular extension, i.e. the diameter of the cold region, a. Clearly, 
a is smaller than 26 l , which is the angle subtended by the full LTB region. The shape of the temperature 
profile (hence a) is determined by the shape of the Void density profile. As for the amplitude A, we recall 
from Q that: 

A « 0.5 Si [LHof (l - ^) = 0.5 ^ . (6) 
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So, A turns out to depend on the radius L and the density contrast Sq, with a dependence on D = L/ tan^L 
which is mild unless the Void is close to the LSS. Note that the expression in eq. © fully agrees with the 
one obtained in (sf, except for our pre- factor which is three times larger, probably due to the different shape 
of the profile (see also appendix We fix the numerical values of A and a phenomenologically. Clearly, 
this leaves a degeneracy in the choice of the physical parameters, since we have three of them (D, L and 
(5o) and only two observational constraints [A and a). In order to fix a range for the numerical values of 
the amplitude A, we proceed as in Q relying on the values given by Q, as follows: for the temperature at 
the centre we use the range T = —(190 ± 80)/iii', which means ^ = (7 ± 3) x 10"'"'; for the angular size a 
of the cold region, wc choose the particular but representative values 6°, 10° and 18° (which correspond 
respectively to Ol = 7°, 11°, 20.5°). 



B. Lensing Temperature Profile 

In this paper, our goal is to compute the contribution AT^^'' /T to the temperature fluctuation, due to 
the lensing of pre-existing primordial temperature fluctuations. The temperature fluctuation due to lensing 
is usually computed in a gradient expansion [ill . , where the total fluctuation is given by 

AT AT(-P) AT(-P) Ar(^) 

— (n) ^ ^^(n) + d,—-{n)d^Q{n) + d,d,^^{h)d^ e{h)d^e{n) + ... . (7) 

We call the correction terms AT'^^^ '^^^ /T and AT^^^ /T, where the superscripts denote the perturbative 
orders. In order to compute this, we need 9, the so-called Lensing potential, which is related to the 



gravitational potential $ by the following line integral 10|, lUl, ll3| : 



V^e = -2 dT^^ — ^V_L$, (8) 

Jtlss '''LSS 

where tq and tlss denote respectively the conformal time at the observer and at the LSS; V_l stands for 
a gradient in the direction transverse to the line of sight. Now, given the gravitational potential, we can 
numerically compute the Lensing potential. The way we proceed is that we approximate the trajectory 
as a straight line and we compute the integral, as a function of the observational angle 6 (in spherical 
coordinates, having fixed the z axis towards the centre of the Void). The gravitational potential <& is zero 
outside the LTB patch and positive at the centre (since we have a Void) . Therefore its transverse gradient 
is negative, which means that the displacement vector is positive (and vanishes at = and 9 ~ 9l). 

Then, we integrate W±Q{9) along the 6 angle and we get 0(0), which turns out to be negative. 

Now, the amplitude of the quantity Vj^0 scales linearly with Sq (since the effect is computed in the linear 
theory) and quadratically with the scale of the structure L: Vj_0 (x Sq {LHq)^ . Removing the gradient, it 
turns out that 9 scales as L^. The final equation for Q can be written as^: 

e{9) = eop{0) , 00 « ^ I'Jol (ii^o)'^ , (9) 

where numerically p(9) is given in fig. [TJ 



^ In the published version of this paper we neglected the dependence on D coming from eq. ||8{. The inclusion of this 
dependence actually goes in the direction of strengthening the results as compared to the published version. 




For reader's convenience, we also give here an approximate polynomial interpolation of the profile: 



p(0) = -l-O.257(^—j +4.076 (^—j - 2.848 (^—j , 9 < 9l ■ (10) 

Clearly p{9) vanishes for 9 > 9^. 

Similarly to the RS temperature fluctuation, also the Lensing fluctuation, AT(n)'^^^/r, is described by 
two parameters, namely those that characterize the Lensing potential &{9): its amplitude Qq at the centre 
of the Void and its angular extension, determined hy 9l- The amplitude Qo can be determined by choosing 
specific values for L and Sq in eq. Note, however, that the dependence of 6o on L and Sq is different 
from the one of A: so it is possible for two Voids to have the same RS effect but different Lensing effects. 
Keeping fixed A as indicated by the observations (A = (7 ± 3) x 10~^) and recalling that D = L/tanS^, 
we can calculate i5o from eq.® by choosing L and 6^: 




2A 



2 tune L 



(LHo)-'^' . (11) 



as shown in the left plot of figl5]for a = 18°, 10° and 6° (from bottom to top). Note that since \So\ < 1, 
one finds a lower bound on L, with a very mild dependence on a, given by: 

i>i„„,.„M^»a5o;-)^. (12) 

This minimum value of L corresponds to the configuration in which the Void is closest to us. 
Wc can put ea. (|TT|) into eq.©, obtaining 

eo^f^^^V" . (13) 



1 



2 tan Si 



as shown in the right plot of fig. [51 This figure shows that the dependence of Qq upon a turns out to be 
not too strong. For instance, the minimum value 8o = 3 x 10^** can be achieved with L « 200 Mpc/h 
for (7 = 6°, while the value Go = 10^^ can be achieved with L = 200 - 500 Mpc/h for a = 18° and 
L = 500 - 800 Mpc/h for a = 10°. 

However, given a fixed value of 9^, both do and Qq blow up at a certain L. This maximum value of L 
corresponds to the configuration in which the Void is very close to the LSS, without touching it. This fact 



6 




L [Mpc/h] L [Mpc/h] 



FIG. 2: Left: \Sq\ as a function of L and for a = 18°, 10° and 6°, from bottom to top. Right; Plot of Qo as a 
function of L. The curves corresponds, from left to right, to cr = 6°, 10°, 18°. In both plots, the shaded regions are 
obtained by varying A in the range (7 ± 3) x 10^°. 



can be understood via the following argument. In a matter dominated Universe the angular size diameter 
distance of the LSS is Da^ss ~ '^/{^lssHo) « 5.45 Mpc//i, where zlss = 1100. Hence, the maximum 
value of L today is 

Lmax ~ ZlSS Da^ss : (14) 

where Ol « L/D. For a = 18° {Ol = 20.5°), a = 10° {Ol = 11°) and tr = 6° {Ol = 7°) this corresponds 
respectively to L^ax about 2100 Mpc//i, 1145 Mpc//i, and 730Mpc//i. These values are in agreement with 
figEl Close to the LSS, the amplitude A is suppressed by the factor 1~DHq/2 of eq.®. Since by definition 
< 1, there are values of L such that the amplitude A does not reach the value 7 x 10~^. For this reason 
there is a maximal value for Gg, as shown in fig. [5] 

On the other hand, when the Void is small, i.e. it is in the position closest to us, reaches its minimum 
allowed value, which is about 3 x 10~^. Clearly, we find this minimum because we are imposing A to be 
in the range suggested by present Cold Spot observations, A = (7 ± 3) x 10~^: had we imposed a smaller 
value of A, we would have obtained a smaller minimum value for Op. 

Note also that assuming a different shape for the profile (see eq. (O) could lead to a different pre- factor 
in eq. dH): instead of 0.5 we could have a pre-factor ranging from 0.05 to 1.8. This would just translate in 
a proper order-unity rescaling of eq. (jl3p and fig. (3] 



C. Decomposition in spherical harmonics 

Given the temperature anisotropy AT'^*^(fi)/r and the Lensing profile 6(n), we will need their spherical 
harmonic decompositions, defined respectively as: 

4L ^Jdn ^Iplll r;^(n) , bi^^J dh e(n) y/„(n) . (15) 

Since the RS temperature anisotropy and the Lensing profile are axially symmetric and since we have 
chosen the z axis to point towards the centre of the Void, the only non-vanishing and bgm are those 

with m ~ and which, in addition, are real. In the left plot of fig. [3] we show the ratio a^^"^^/^ as a function 
of the multipole £. From bottom to top, the curves correspond to a Cold Spot with diameter a = 18°, 10° 
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and 6° respectively (for more details see p). In the right plot of fig. [3] we show the corresponding bio/Qo 
coefficients of the Lensing potential 0(n). 




70 i 



FIG. 3: Plot of a^Q /A and beo/Qo, as a function of the multipole £ for a = 18°, 10°, 6° (from bottom to top). 



Given the b^ coefficients in (|15p . we may compute the first order a[^^^^ coefficients for the Lensing 
temperature profile AT^^'^^^/T, via the expression [lo| : 

AL)ii)_ V r— 'Of ^yn+m' i'i^' + !)-£{£ +!)+£"{£" + 1) jp). 

«£m ~ 2^ e ('£"(.- ^) 2 a^,_^,bi„o , (15) 

e',m',e" 

where we have also used the fact that ben m" oc Sr,i"o, and that they are real numbers. We have introduced 
the Gaunt integrals, which arc given in terms of the Wigner 3-j symbols [lo| as follows: 



r.m,m,ms - , / (^^1 + + 1)(2^3 + 1) f h h 4 \ f h h 4\ 

^2 «3 "V 4^ 1^0 O^Vmi 7712 ms) ' ^'''> 

The Gaunt integrals are non-zero only if the sum of the upper indices is zero, so we can set m' = m in 
ea.([T6| obtaining: 



(L)(i) X- r^-™™o i'ie' + i)-eii + i) + rii" + i) (p). 

'^m - ^ f 5 af'_„O£"0 • (,tOj 



We recall that (a)„, ) = 0, where the brackets stand for a statistical average over an ensemble of possible 
realizations of the Universe - or, equivalently, an average over many distant uncorrelated observers. Given 
ea. ([T5|) . it turns out that the same applies to the lensing coefficients, namely (a^^^*'^'') = 0. 
The expression for the second order coefficients a^^^ ^^'^ has been derived 6.17. in 14|, eq. (46). 



III. TWO-POINT FUNCTIONS 



Given a temperature profile with its aim coefficients one can compute the associated two-point correlation 
functions. In general, given a single set of aim coefficients, the two-point correlation function (power 
spectrum) is defined via the C^'s coefficients as 

y \aiml_ ^ 
^ 2£ + 1 ^ ' 
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Note that this definition ensures that the CiS do not depend on the choice of the coordinate system on 
the sphere. Therefore, as in Q, we may keep our z axis ahgned with the centre of the Void. In our case 
aim = 0/irn ^ 'Hm ^-nd, in order to face with the experimentahy observed values for Ci, we have to 

estimate the theoretical prediction for (C^). 

For a primordial Gaussian signal the two-point correlation functions are given by: 

(4fi«lfi:>='5....'5„,„,(c(f) , (20) 

where the (C) ) are predicted by some mechanism {e.g. inflation) that can generate primordial Gaussian 
fluctuations. In Q we already computed the contribution to the power spectrum given by the RS effect, 
Here we consider the Lensing term, which leads to additional contributions of the form 
(PL contribution) and (LL contribution). The PL contribution actually contains 

two terms: one is first order in the Lensing potential, (a^'^-'a'^-' '^•'*), and one is second order (a^'^^a^^^ (^)*). 
It turns out that the latter term is of the same order as the contribution and with opposite 

sign, leading to a cancellation^, as happens in Lensing from usual primordial Gaussian profiles 13 1. 



Given the expression in eg. ljlSp . we may compute the first order contribution to the two-point correlation 
function due to the primordial and Lensing temperature fluctuations, as follows: 

, [P] (L)(l)*. _ r-r»2m20 + 1) ~ ^2(^2 + 1) + + 1) , .yn^iiP) JP)*\h 

f ,1" 

I" 

The diagonal term leads to the following contribution to the C^'s, 

t I (P) 

which turns out to vanish, as one can see substituting the explicit expression for the Gaunt integral and 
using a well known property of the Wigncr 3-j symbols: 

<c<-x..,Hcr>i:;^(^ I :)nr+i,v.i;(-ir(_L /„ :)^o, ,23, 

^ V ' 

The non-diagonal terms are quite interesting: in fact ([2T|) leads to a correlation between different £'s. The 
correlations arc small, but they are present also at high ^'s. In fact any ^\ and ^2 will be correlated as long 
as \i\ — ^2! ^ 60, because the Gaunt integrals in ([?!]) are nonzero if \l2 — l\\ < i" and the coefficients 
are non-zero for ^ < 60 (for the values of a relevant for the Cold Spot, see fig. [3]). This signal could be of 
interest when looking at high-resolution experiments such as Planck, because it could reveal the existence 
of the large Void through the Lensing effect. An interesting formalism to analyze these correlations in a 
rotational invariant way would be to introduce the so-called "Bipolar power spectrum" [l^, which is a 



combination of the non-diagonal a^m's, which generalizes the C^'s. The Bipolar power spectrum is zero 
for Gaussian fields, but it has a contribution due to cosmic variance. It would be interesting to apply this 
analysis to our case and we postpone this to future work. 



^ We thank A. Challinor for pointing out this fact. 
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The pure Lensing contribution to the power spectrum, defined as 



, (L)(l) (L)(l)*v 
\"^m "'irn I 



2£+l 

m=-l 

upon substitution with the expression in eq. (jl8p . turns out to be 



(24) 



E 



2£i + l 
47r 



£i(£i + l) - £(£+!)+ ^2(^2 + 1)' 



o)'(^^^o)' ■ (25) 



^ ^1 I 




Note that this expression agrees with the last term of cq. (62) of IJ], using the fact that the power 
spectrum of the Lensing potential in our case is 6^g/(2^ + !)■ In addition to this term we have to add the 



term (0^^ a*-^-* ''^''), which is also given in eq. (62) and (63) of 14 1, and which is equal to 



where 



which can be easily computed and gives the result 

Therefore the final result for the correction ACe to the power spectrum is given by [l^ [l^ '■ 



(26) 



(27) 



(28) 



The first term is a multiplicative correction to the primordial power spectrum and the constant i? can be 
interpreted as the average of the square of the deflection angle [13| . The pcrturbative treatment of Lensing 

+ l)i? < 1. This means 



is justified as long as this correction is small with respect to {Cg ), that is 



that the approximation is good until some 



which is given in our case, by: 
^ _ 2.5/60. 



This is a good approximation for £max ^ 2000, provided that 9o S 10 ^- In order to go to higher 
resolutions one should include the Lensing effect in a non-perturbative way, see e.g. |13|. 

The second term in eq. (|28|) is a convolution: it can transfer power from large to small scales and in 



general can smooth some of the primordial peak structure. The two terms arc similar in magnitude but 
with opposite signs, and they cancel out with a precision of about lO"'^, leaving an oscillatory correction. 
The final result for the ratio ACi/ {C'^^^) is shown in the left plot of fig. 2] as a function of the multipole 
i, assuming a specific value for the amplitude of the Lensing potential, namely the minumum allowed one, 
00 = 3 X lO^"'. Note that for different values of Go one has simply to multiply the results in fig. H] by the 
factor 6^/(9 x lO^^). 

For an experiment whose sensitivity goes up to some (.max, the signal-to- noise ratio associated with the 
pure Lensing contribution is given by: 



LL 



E 



AC? 



(P) 



where uf is the cosmic variance of the power spectrum: aj ^ 2/(2^+ 1) {C^ 
contain also the experimental noise, but we have just cutoff the sum at the t 



(29) 

^. In general should 
corresponding to the 
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1000 2000 3000 4000 5000 1000 2000 3000 4000 5000 



FIG. 4: Left: Plot of ACt/{C^/'^) as a function of the multipole £ and with the minumum allowed value Go = 
3 X 10^*. Right: Plot of {S/N)ll as a function of the multipole £. The shaded region is obtained by varying A 
in the range (7 ± 3) x 10"''. For the lower curve we have taken L = 200Mpc//i and a = 6°, which corresponds 
- see figH]- to 0o — (S.l^'y) x 10"'*; for the middle curve we have taken L — 400Mpc//i and a — 10°, namely 
So = (7.4l};^) X 10""; for the upper curve we have taken L = 800Mpc//i for a = 18°, namely So = (2±.5) x 10"^ 
Both curves are a good approximation only up to ^max ~ 2.5/0o, while for higher £ there is an C(10%) error, due 
to the fact that the gradient approximation breaks down [l^ . 



experimental sensitivity. The signal-to-noise ratio is displayed in the right plot of fig. U for various Void 
configurations. Note that {S/N)ll scales as 6g. 

Fig. U shows that, if one assumes that the CMB Cold Spot is due to the RS effect of a Void with central 
temperature fluctuation yl = (7± 3) x 10"^, then very high-resolution CMB experiments could detect Qq: 
WMAP could see a signal only in the extreme case if Qq is bigger than about 3 x 10"^ (which corresponds to 
a cr of at least 10° and L > 1000Mpc//i, see fig. [2]), while Planck could go down to 10"^, which corresponds 
to L > 500Mpc/ft. for a = 10°. In order to be able to detect smaller Voids, one should look at the Cold Spot 



region with a higher resolution experiment (such as the planned Atacama Cosmology Telescope 15| or the 
South Pole Telescope 0): e.g. going up to ^ ~ 5000 seems necessary to detect Voids with L ~ 300Mpc//i 
and (7 = 6°. 

Moreover note that, clearly, if there is such an increase in the power spectrum due to a single Void, 
this would be localized in the Cold Spot region: therefore, masking this area the effect should be absent. 
It is worth also mentioning that there could be other Voids of similar, or smaller, size in the sky, (for 
example (l?} points out out that there might be other cold and hot spots of similar size in the Southern 
hemisphere). This would lead to similar contributions to be added to the power spectrum, with different 
amplitudes. 

An intriguing consequence of our analysis is that, if the Void is large enough, its Lensing signal could 



18, hJI in the 



possibly be linked to the recent claimed detection of a hemispherical power asymmetry 
WMAP data, which extends also to high multipoles {£ ~ 2 — 600). While the \ow-£ power spectrum 
asymmetry could be linked to the existence of large structures (due to the RS effect, as discussed in Q), 
the fact that the signal is not limited to low multipoles could be explained by the lensing effect of such 
objects. Indeed, the direction of the center of the hemisphere with the largest power happens to be close 
to the location of the Cold Spot. The amplitude of the claimed effect [l^ is of order 6Ce/Ci ^ 10"^, which 
would be compatible with our effect (see fig. H]) ii &o is of order of a few times 10"^, which corresponds to 
L ^ 1000Mpc//i. Therefore, it is observationally important to determine whether there is an increase in 
small-scale power, with the shape we predict in fig. 21 which could be linked to a huge Void, or possibly to 
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several Voids and structures on large scales in the sky. 

We may, finally, mention that the recent paper Q has proposed to estimate the amplitude of the lensing 
effect by looking directly at AT^-^^/T, instead of looking at a quadratic estimator such as (C^'^'^-'). This 
leads to the consequence that the signal-to-noise ratio is linear in the (c]^^'')'s, instead of containing its 
squares, as happens for our eq. ((29|) . The maximal value of the lensing potential amplitude of Q coincides 
with our 60 for L = 350Mpc//i and = 0.3, namely 60 = (7 ± 4) x 10"*. Note however that [i| is 
considering a different geometry: a non-compensated cylinder of radius 150 Mpc and transverse length 
200 Mpc. The signal Si in cq. (13) of Q is similar in structure to one of our terms, the (C^^^"*) in eq. (P5|) : 
in fact they turn out to have a very similar shape, although our (C^'^'^'')'s are smaller by a factor of order 
unity. However, the existence of an additional correction of the type (a^^^ a^^-"^^-') may affect, as in our 
case, the conclusions of since it might lead to a cancellation with a precision of order 10"'^. 



IV. BISPECTRUM 



After having computed the coefficients and the two-point correlation functions we may estimate the 
impact of a large Void on the bispectrum coefficients. In [6| we have considered the bispcctrum due to the 
RS effect only, while in the present paper we focus on effects due to lensing. 

The basic quantities arc the BJ^'-'T^^P^ coefficients, defined as: 



fi £2 is — "^1"'! ^^2^2 a^sma ■ I-JLIJ 

These quantities are coordinate-dependent, so we define (in analogy with the C^'s coefficients), the angularly 
averaged bispectrum 

Be^e^e, = Y f ^1 ^2 ^3 \ 5-i'"2m3 _ (3^^ 
tit2t3 \ mi 7712 ma y ^ ^ 

One can show, using the properties of the Wigner 3-j symbols [lo| . that this combination does not depend 
on the chosen z-axis, so these are the quantities that we can use to make predictions. For convenience we 
keep our z-axis along the direction of the centre of the Void. Then we perform the statistical averaging 
of the primordial perturbations. As we have discussed in Q, using eqs. (PH)) and (|15p we get a sum of 27 
terms. However, we want to give a prediction, i.e. compute statistical averages (Bi-^i^i^), and several terms 
have zero average. As we have said, a crucial assumption that wc make is that the coefficients o^q are not 
stochastic quantities, which means that the location of the Void in the sky is not correlated at all with the 
primordial temperature fluctuations coming from inflation. Note that this is a conservative assumption: 
if there is some correlation, and barring cancellations, the non-gaussianity could be much more important 
since some terms would be non-zero. Under this assumption, we have shown in Q, that the two leading 
terms are 



,m2 ^ ^ 



and 



?ni ,m2 ,m3 

The first one has been studied in detail in Ifil- We now consider the second. 
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A. Signal-to-Noise ratio 



In this subsection wc compute the contribution to (Bf-^i^t^) due to the couphng between the Primordial- 
Lensing and RS effects, as fohows: 

(^rrr)™^' = /(1,2,3) + /(2,1,3) + /(1,3,2) + /(3,1,2) + /(2,3,1) + /(3,2,1) , 



(34) 



Substituting the exphcit expressions for the Lensing coefBcients, ea. (jl8p . we get: 

e',e" 

where we have used the fact that (a^fii^ a^r_\*„2 ) = (C^'f ■')'5^i^"5mi-Tn2- 
The angularly averaged bispectrum is given by: 



\ «i«2*!3 / V TOi TO2 J ^i^2<:3 / ' 

and can be written as 

(■^wT') = ^(1, 2, 3) + F(3, 1,2) + Fi2, 3, 1) + (-1)^1+^2+^3 (^(2^ 1, 3) + i^(l, 3, 2) + ^^(3, 2, 1)) , 

(37) 

Substituting the expression for the Gaunt integral in (|35p we find 



(P), (PS)Y- /(2£j + l)(2£. + l)(2r + 1) e"\ k{k + 1) - ^j-(£j- + 1) + + 1) 



47r 10 



so that, finally. 



FU . k) /r(^)\.(«^) / (2^, + i)(2£. + i) fe, 4A £.(£, + i)-4-(4- + i) + 4(4 + i) , 

nhJ,k)-{C,^ }a,^, y 4^(24 + 1) 1^0 oj ^ " ^'^^> 

It is customary to define a reduced bispectrum 6fif2^3 the following: 



_ (24 + l)(2£2 + l)(2£3 + l) / 4 4 4 1 , ,4.^ 

We may then easily compute, using the {b^ut^^^) coefBcients: 



{f>ui ) - {^m ) (2^^1)3/2 l^OOoj "2^ ^ ' ' 2e+l ' '° 
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FIG. 5: The (6^^^ ) coefficients - with ^ even - due to Primordial-Lensing-RS for angular diameter cr — 18°, 10°, 6° 
from left to right. We assumed a Void with Oo = 3 x 10^*. The solid line shows the primordial contribution with 
/nl = 1. 



Note that {h'^ui'^^^) = if ^ is odd. The result for the {b^^/^^^^} with i even is shown in fig. O taking 
6o = 3 X 10""* and, from left to right, a = 18°, 10°, 6°. The {h^^t^^'^) coefficients look comparable to the 
primordial non-Gaussianity with /at^ of order unity (the solid line is obtained with /at^ = 1) and they are 
non-zero for low t. 

As can be seen by inspecting (|37l) and (p9)) . {B^^^.^ff'^) is non-zero provided that one of the i indices is 
low, say £ < 60 for cr > 6° (see figEl). This means that, for an high resolution experiment as WMAP or 
Planck, which goes up to very high there are many non-zero coefficients, with Ik ^ 60 and arbitrary li 
and ^j. The only constraint on £i and i.j is given by the triangular condition of the Wigner 3-j symbols 
(which vanish unless \ti ~ lj \ < < + for any £j and Ik). 

Wc may indeed compute the Signal-to-Noise {S/N) ratio for the PLRS efi'ect. For a signal labeled by i, 
this is defined as 



10|: 



{S/N)., = ^ , Fu= E (41) 

^JPr:^ 2<e^<e2<es<i^^^ '^^i^2«3 

Here crfj^^^g is the cosmic variance of the bispectrum 

^?if2f3 ^ {Ctl){Cl2){Cl3)^llt.2t.3 , (42) 

where A^j^^^, = 1,2, or 6 respectively if all ^'s are different, if only two of them are equal or if they are 
all equal. The {Ce)^s are the sum of the CMB power spectrum plus the power spectrum of the noise of the 
detector. In general at some iraax thc uoisc becomes dominant, while below it (C^) is dominated by the 
primordial noise: {Ci) ~ {C} ). 

The result for {S/N)plb.s is shown in fig. [6l assuming that the detector noise is negligible. Note that 
this ratio scales as the product AQq and is nearly independent on cr. Indeed, even though the coefficients 
a[|f'^^ and 6;o display a dependence on the angle cr (see fig. [3]), we have found that the dependence of 
{S/N)pLjis on a is very weak for 6° < ct < 18°. Keeping A fixed - the shaded region corresponds to letting 
A vary in the range (7 ± 3) x 10^^ - the Tensing signal is clearly bigger the larger Op (hence L) is. From 
left to right, the three panels in fig.[S]show {S/N)plbs in the case that L = 800,400,200 Mpc//i (we have 
chosen for illustration various values of cr, corresponding to Gq = (2 ± .5) x 10"'^, 9o = (7.4j^i[ g) x 10~^, 
00 = (3.ll;7) X 10^^. One can see that this signal is qualitatively different from the signal induced by the 
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L = 800 Mpc//i L = 400 Mpc//i L = 200 Mpc//i 




FIG. 6: The Signal-to-Noise ratio for the Primordial-Lensing-RS couphng. The shaded regions are obtained by 
varying A — {7 ± 3) x 10"''. For the right plot we have taken L = 200Mpc//i and a = 6°, which corresponds 
- see figO - to &o = (3.ll'7) x 10"'*; for the middle plot we have taken L = 400Mpc//t and a = 10°, namely 
Bo = {7-4:l]i) X 10""*; for the left plot we have taken L = 800Mpc//i for a = 18°, namely Oo = (2 ± .5) x 10"^ 



pure RS contribution (see fig. 6 of ref. [Gf) which becomes constant for imax ^ 50. In fact, now the effect 
is tiny for small imax, but it is sizable when £max is large. As a result, the total (S/N) plb,s increases with 
the sensitivity of the experiment, the physical reason being that any i is lensed by the Void, also the small 
scale ones. 

As it can be seen from the plot, the ratio {S/N)plb.s summed over about 1900 multipoles (the maximal 
resolution for the Planck satellite) is significantly larger than unity, roughly if i > 300 Mpc//i. We stress 
therefore that this signature can uniquely discriminate the presence of a Void: in fact, if the Cold Spot is 
just a statistical fluke and if there is no Void (or a similar object) in the line of sight, the Icnsing effect 
would be absent. It is also clear from the plot that going to slightly higher resolutions (as can be obtained 
by 15[ or 16|) will probe the entire parameter space. 



B. Contamination of Jml measurements 



We want here to find out what is the impact of having a hypothetical huge Void in the line of sight 
on the measurement of the primordial non-gaussianity parameter Jml- It is customary to analyze the 
non-Gaussian signal via the quantity Jnl- Its definition is given [lO| for the so-called non-gaussianity of 
local type, parameterizing the primordial curvature perturbations 0(a;) as follows 

d^{x)=<j>L{x) + fNL{4>l{x)~{4>l{x))) , (43) 

where (/>l(x) is the linear Gaussian part of the perturbation. Generically, given a physical model {e.g. slow- 
roll inflation), fj^L is a function of the momenta, but it is usually assumed in quantitative data analyses 
to be just a number. Single field minimally coupled slow-roll inflationary models predict generically very 
small values, /at^ = O(O.l) l3, 2^, 21 1, while other models may predict larger values [lol |. 



As discussed in [6[, the RS effect leads to a large correction to (Bi-^i^i^) for 2 < £ < 60. The Lensing 
effect, instead, is much smaller at low-£, but it can contaminate the signal at large £, since it couples the 
low-^ of the RS effect with the high £ of the primordial signal. Wc now estimate the amount of such a 
contamination. 

The primordial bispectrum coefficients can be written as [lo| : 

^Ill2t3 - ^i^i^i^ , (44) 
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where the have a specific form in terms of the primordial spectrum and the radiation transfer 

function 



io|. 



The impact on /jv l due to the PLRS effect can be computed estimating the following ratio 10|, l22[ : 



A/ir/^^)(€™a.) = . (45) 

2<£i<£2<^3<f™a. <T^ , ' 
'1'2*3 

In general this is computationally very demanding, since it would require to compute the primordial i?^^"^^ 
up to very high multipoles. However, the effect due to the Void is dominated by squeezed triangles: since 
one side is given by the RS effect which goes up at most to ^ 80, in order to compute the cross-correlation 
in the numerator of ea. (|45p we only need the squeezed configurations also for the primordial effect. We may 
give a rough estimate, using the so called flat-sky approximation, in which one of the ^'s is much smaller 
than the other two. In this case we can write down a very simple expression for the coefficients, 
following [23j . As discussed in Appendix [B] the result of this approximation turns out to be: 

6?;;";^«-i2(c(f))(c(f) . (46) 

The denominator of eq. (j45p is instead a given quantity for a given experiment. In fact, the experimental 
noise in (jf^^^^^ at some high Imax (dependent on the experiment) becomes so large that the multipoles 
^ > ^max do not contribute to the sum. The denominator in ea. ps)) is roughly equal to 10^^ x i'^^.x^ ^ 



shown in [10 1 . 



Performing the sum we get a result which is of order of 



Using the maximal value for A and Op, we can go up at most to 10^"^. The result has some imax dependence, 
but it is always of the order of 10~^, even for Imax ^ 2000. 



A \ / e 



V. CONCLUSIONS 



Motivated by the so-called Cold Spot in the WMAP data, we have studied in this paper the impact on 
CMB two and three point correlation functions of the presence of an anomalously large Void along the 
line of sight, whose existence could be at the origin of the Cold Spot. In particular, we have extended the 
analysis performed in Q by including the Lensing effect on the CMB primordial photons. 

For the power spectrum, the Primordial-Lensing coupling vanishes exactly. It leads, however, to non-zero 
correlations between different ^' s in the off-diagonal two-point correlation functions. The Lensing-Lensing 
coupling, instead, gives a non-zero signal in the power spectrum, which turns out to have a Signal-to-Noise 
ratio of order unity already for WMAP sensitivity only if the Void radius L is extremely large: bigger than 
about 1 Gpc//i (since the amplitude of the Lensing effect increases with L). The Signal-to-Noise ratio 
turns out to be larger than unity for Planck sensitivity for L > 500 Mpc//i. 

In addition, we pointed out an intriguing consequence of our analysis: it could be possible that the 
existence of one (or several) large scale Void(s) in the sky appear not only as anomalies at low ^, but also 
at high I, through the lensing effect. In the case in which the Void is very large (i.e. L ^ 1 Gpc//i), this 
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could be linked to the claimed detection of a hemispherical power asymmetry in the WMAP data [18| . 
which extends also to high multipoles = 2 — 600) 19[. 

For the bispectrum, there is a non-zero coupling between the Primordial fluctuations, the Rees-Sciama 
effect (present at low t) and the Lensing effect (present for any €). This would provide a distinctive 
signature in favor of a Large Void. In fact, its Signal-to-Noise exceeds unity for the Planck experiment if 
the radius of the Void is at least about 300 Mpc/ft-. Higher resolution experiments will be able to detect 
this non-gaussian signal also for smaller Void radii, covering the entire parameter space. 

Finally, we have studied the impact of such a structure on the determination of the primordial non- 
gaussianity, finding that the contamination is negligibly small. 
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Note added: In the published version of this paper we neglected the dependence of the Lensing Potential 
amplitude on the comoving distance of the Void, see eq. The inclusion of this dependence actually 

goes in the direction of strengthening the results discussed in the published version. 



Appendix A: Inclusion of a cosmological constant 



In this appendix wc show how our results change by including the effect of a cosmological constant, 
which has the qualitative effect of increasing the RS effect because a net redshift is already present at the 
linear level in 01 • We remind the reader that in this case the effect usually goes under the name of 
Integrated Sachs- Wolfe effect (ISW). Therefore this additional ISW effect is especially relevant when |Jo| 
is much smaller than unity, which is the case if the Void radius is very large, but still not close to touch 
the Last Scattering Surface (as can be seen in fig. [5]). 

For this estimate we use the calculations performed by Q , with the caveat that they are performed in a 
different setup, with a different density profile, namely in the approximation in which the Void is uniformly 
underdense and it is compensated by an overdense thin shell. In the limit of zero cosmological constant 
our expression eq. ^ fully agrees with the amplitude derived in 0|, except for our pre- factor which is 
three times larger, probably due to the different shape of the profile. We use the expression, multiplied 
by a factor of 3, given in [5[ for the amplitude A of the (RS)-I-(ISW) effect in order to derive a relation 
between the density contrast at the centre of the Void and the radius of the Void (in such a way that 
A = (7 ± 3) X 10^^). This relation is shown in the left plot of fig. [71 choosing the matter density fraction 
Hm = 0.27. This can be directly compared with the left plot of fig. [21 for which ~ 1. Then, we 
can plug the new relation for 5q into our eq. (O, which allows us to reconstruct the value of L given a 
measurement of the Lensing amplitude Qq. We show the result for VIm = 0.27 in the right plot of fig. [71 
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500 1000 1500 2000 500 1000 1500 2000 



L [Mpc/h] L [Mpc/h] 

FIG. 7: Left: \So\ as a function of L, with fijv/ = 0.27 and for a — 18°, 10° and 6°, from bottom to top. Riglit: 
Plot of 00 as a function of L, with Qm = 0.27. The curves corresponds, from top to bottom, to cr = 6°, 10°, 18°. 
In both plots, the shaded regions are obtained by varying A in the range (7 ± 3) x 10"''. 



which can then be compared with the right plot of fig. [2] for which Hm = 1- 

For illustration, in fig. [8] we also plot the ratio of 6o obtained for VIm = 1 over the one with VIm = 0.27. 
As we have anticipated, the ratio is larger than 1 in the regions of parameters where Sq is much smaller 
than unity. Anyhow, we stress again that this ratio is only a rough guide, given the fact that 0] is using a 
different profile. 




L [Mpc/h] 

FIG. 8: Ratio r defined as Oo with Oa/ ~ 1 over Go with Qm = 0.27. The curves correspond to cr = 6°, 10° and 
18°, from bottom to top, and to A = 7 x 10"^. 



Appendix B: Flat-sky approximation 

We re-derive here the approximation for the bispectrum , following [2^ . As a starting point we take 
the expression: 

C^^a = J r'dr[f3e, {r)P,, {r)ae, (r) + (r)a,, {r)f3e, (r) + a,, (r)/3,, {r)f3,, (r)] , (Bl) 

where 

Mr) ^ ^ f k^dkP{k)ji{kr)Ai(k) , (B2) 
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and 



ai{r) - - / edkn{kr)Ai{k) . (B3) 



TT 



Here ji{kr) arc the spherical Besscl functions, P{k) is the primordial power spectrum and A.i{k) is the 
radiation transfer function (which goes as A;(fc) w ^jiikr^)/?! for large scales). Defining tq as the present 
day value of conformal time, tr as the value at decoupling and = tq — t^j as the comoving distance to 
the surface of last scattering, the region of integration for r is over the sound horizon (from tq to rg — 2tr). 
We can rewrite (IBll) as: 



C^™3 = ^/^i / fc?dfcifc2^^fc2fc|dfc3A,,(fci)A,,(fc2)A,3(fc3)C,,,,,3(fci,fc2,fc3)[^(fci)/'(fc2) + c2/c.] , (B4) 

where we define Cg^ir^g^{ki,k2,k^) = j r'^drji^{kir)ii^{k2r)ji^(k^r). In the limit of collapsed triangles, 
t-i ^ ^2j^3, we can approximate 



Ceie^e^{ki,k2,k3) ji^ikiro) J r'^drji^{k2r)ji^{k3r) , (B5) 
and use the following property of the Besscl functions 

'r2drj,,(fc2r)j,,(fc3r)^|^fc-^ . (B6) 
Substituting this result into (jB4|) we find 

^ZZ^^-^Inl J fc?dfcifc|dfc2if,(fcirz,)A,,(fci)A,,(fc2)A,,(fc2)[P(fci)P(fc2)+c2/c.] . (B7) 
This can be evaluated as 

bi,i,e,^-UfML{d,P}{C^,p} , (B8) 
where we have used the fact that )c[^^) is given by: 

{cP) = / dke\Ai{k)\P{k) (B9) 



{2nf 

and where we have neglected the last term in the cyclic permutations, which is proportional to P(fc2)^, 
since P{ki) > ^(^2). 
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